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Abstract 

We present a dynamical model for the price evolution of financial assets. The model is based in 
a two level structure. In the first stage one finds an agent-based model that describes the present 
state of the investors' beliefs, perspectives or strategies. The dynamics is inspired by a model for 
describing predator-prey population evolution: agents change their mind through self- or mutual 
interaction, and the decision is adopted on a random basis, with no direct influence of the price 
itself. One of the most appealing properties of such a system is the presence of large oscillations in 
the number of agents sharing the same perspective, what may be linked with the existence of bullish 
and bearish periods in financial markets. In the second stage one has the pricing mechanism, which 
will be driven by the relative population in the different investors' groups. The price equation will 
depend on the specific nature of the species, and thus it may change from one market to the other: 
we will firstly present a simple model of excess demand, and subsequently consider a more elaborate 
liquidity model. The outcomes of both models are analysed and compared. 
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I. INTRODUCTION 



Financial models based on interacting agents possess a large tradition in the economic 
literature [1'] — one of the first references in which the evolution of a market is related with 
the activity of individual investors dates back to 1974 — , but they have gained relevance in 
the physics literature in relatively recent years . The complete list of such models is so 
extensive and their properties so diverse that we can merely sketch here the recurrent traits 



shared by most of the models, and address the reader to the references cited in 

The pioneering work of 1974 ^1 already contains one of the more ubiquitous ingredients 
in the subsequent agent-based models: heterogeneity. Agents are assumed to be heteroge- 
neous to some extend, and therefore they can be aggregated into one out of a finite set of 
categories. Since the minimum number of different categories is two, and simplicity is often 
a plus, investors are usually arranged into two (competing) groups. The terms used to name 
them and their defining properties are uneven across the literature — Chartists and Funda- 

n n n 

mentalists j^. Trend-followers and Contrarians [5], Speculators and Producers p. Imitators 
and Optimizers {7] — but the ideas lying behind are similar, and can be well represented by 
Chartists and Fundamentalists. Chartists are (sometimes adaptative) agents whose inver- 
sion strategy is based in the belief that past information may contain clues about the future 
evolution of the security, and therefore that they can infer future prices. Fundamentalists 
are in essence agents that think that they can deduce the present value of a firm on the basis 
of the information presently available, such as dividend payments or earning rates. Funda- 
mentalists operate in a rather predictable way since they expect that the market tends to 
correct any observed deviation between fundamental and market prices: they sell overpriced 
securities and buy underpriced ones. The picture is not so simple for Chartist-like investors 
since at the end they deploy rule-of-thumb strategies, like the Moving Average Convergence 
Divergence (MACD) indicator or the Relative Strength Index (RSI), two technical analysis 
tools broadly used by actual financial practitioners. Therefore, the list of available strategies 
in agent-based models may be so large that, in the most extreme situation, strategies may 
differ for any pair of investors in the market, as in some instances of the Minority Game 
market model j^, 8|. In fact, any single agent may combine technical trading rules with 
fundamental ones, or decide among them, what makes evolve the profile of the investors. 
There is no doubt that this diversity adds more heterogeneity into the model. 



2 



Another general trait of present models is that agents are boundedly rational j^: they 
decide their actions in the next time step on the basis of a limited and possibly incomplete 
amount of information. They ignore the beliefs of the rest of investors and usually cannot 
evaluate the consequences of their own decisions. Under these circumstances selfish agents 
try to maximize a pay-off or utility function, a measure of their individual success. 

The final ingredient is the pricing mechanism. The usual paradigm when agent activity 
does not explicitly settle the price of the asset is to define a differential equation or a finite 
difference equation that relates the price evolution with the relevant (global) variables of 
the model. Since these variables are affected by the mutual interaction of the investors in 
a complex way, two complementary approaches are generally considered: the behaviour of 
the system is computer simulated and/or the complexity is reduced by considering that the 
number of agents approaches to infinity, the thermodynamic limit. 

As we will shortly show, some of the previous ingredients either are not present in this 
agent-based model, or have been introduced with a different philosophy. The model was 
inspired by a previous article on population dynamics [10], where authors reported presence 
of large oscillations in the species densities due to a finite-size stochastic effect. We export 
this idea into the financial language with the confection of a model that describes general 
aspects of investors' dynamics. The behaviour of the asset price is first obtained after 
assuming a simple model of excess demand, and subsequently, we apply the same approach 



to the modelling of the interaction between limit and market orders in a stock market 



overall result exhibits similarities with prevailing agent-based financial models Uj, 
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The paper is structured in three main Sections. Section [TTl deals with the agent model 
strictly speaking: First we define who are the agents and the three different states in which 
they can be found at every moment, the mechanisms that govern the changes from one state 
to the other, and the transition rates between states. Then we derive a master equation that 
characterises the time evolution of the system, and analyse the stationary solutions of this 
equation in the thermodynamic limit. Finally we find the second-order corrections and show 
their relevance in finite-size models. In Section lXTTl we establish a first connection between the 
agent model and market price changes: we simulate the time evolution of the system under 
representative market conditions, analyse the most relevant traits, and compare them with 
well-known empirical properties of actual financial time series. In Section IIVI we propose a 



3 



second identification for tlie species categories (liquidity providers and liquidity takers), and 
a different price formation procedure is considered. The outcome presents new properties 
that are still consistent with what one may find in practice. This reinforces the potentials of 
the model. Conclusions are drawn in Section IVj and some technical derivations concerning 
correlation functions are left to an Appendix. 



II. AGENT DYNAMICS 



As we have just stated, this Section deals exclusively with the intrinsic features that the 
agent interplay generates. For this aim, as we will see, we do not need a detailed description 
of the agents' inner properties. The most important point to be made here is about the 
motivation and plausibility of the agent-based approach introduced. 

Along this article we will assume that any trader that may ever operate in our finan- 
cial market can be accommodated in two great, well-defined and excluding groups. The 
first and most populated group of investors will constitute what is usually termed as noise 



traders 
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16|. We will assume that each one of these traders acts in a purely random 



fashion, independently of the rest of agents in the market. We are not considering this kind 
of traders as individuals, they merely act as a some sort of thermal bath or noise source, 
what reinforces the foundations of the stochastic character of the dynamics to be introduced. 

The second group of traders is the set of those which we will call qualified investors, but 
the term informed traders [3] would be fit for them as well. As we are showing below, we 
will connect the price evolution with the collective state of these players, so we will consider 
that this group embraces the main actors with bigger influence in financial terms: mutual 
funds, investment banks, or corporations in general. The total amount of such participants 
in a real market is much more moderate [3] , what makes sense to consider a finite-size agent 
model to describe them. We will assume that the instantaneous state of these investors will 
fluctuate within three categories as a consequence of the interaction with the rest of agents. 

Thus, in essence, our agent model will describe how populated are these three states as 
a function of time. The detailed way in which the states ought to be defined will depend on 
the specific financial market under analysis. Different markets may be sensitive to different 
aspects of the agents' activity. We will delay the definitive profiling of the categories until 
the market dynamics is considered, in Section IIIII and Section IIVI and centre our attention 
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now in the species self-interaction, as announced. 
A. The species self-interaction 

Let us consider then a finite set of N fuUy connected interacting agents who, at every 
instant of time t, may be found in one out of three possible states that we will label by 
the letters A, B and E. In a very general sense, that must be further refined from case 
to case, we will assume that an abundance of agents in state A, A^a(^), is related with a 
bear market scenario, that the increasing of population in the B side, A^s(^), leads to a bull 
market scenario, whereas the market is not sensitive to changes in the number of agents of 
type E, NE{t), beyond the fact that = iV^(t) + NB{t) + NE^t) is fixed. We can think for 
example that these three sets categorise those agents that are either willing to sell, to buy, 
or that keep a veiled attitude, but there are several alternative identifications for groups A 
and B — liquidity providers and liquidity takers, or pessimistic and optimistic investors — 
that would equally help in visualising these two groups. In any case, within our formulation 
E will always represent a neutral position. 

Another ingredient of the model is the fact that the mechanism which allows agents 
to change their minds and move from some state to another is based on self- and mutual 
interactions. Note specially that decisions are not affected by the previous history (there is 
no necessity in selling after buying, for instance) what confers the process with the Markov 
property. 

In a population dynamics language, one has two species living in a finite world: As we 
will see afterwards with the description of the rules of engagement, within this model, the 
A^s will play the role of preys and the 5's will be the predators. The E^s, those agents 
without a definite or explicit intention, act as empty space. 

The basic unitary interaction in population problems is the death process, A ^ E and 
B E. Each one of these two processes, and the same applies for the rest of interactions, 
may encompass the aggregate effect of disparate contributors. ^ So, within an offer /demand 
framework we ought to include here all the transactions between the agents and the noise 

^ Note however that our approach for the agent dynamics is mostly phenomenological. The available 
interactions were selected with the aim of capturing general macroscopic features. Therefore the proposed 
and future interpretations of the processes in terms of microscopic events fulfil mainly illustrative purposes. 
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traders, perhaps the most common transaction instance in actual stock markets, but one 
can also take into account the possibility that qualified investors may decide to dechne in 
their previous intentions by their own. In this scheme p and q represent the probability per 
unit of time that a given active agent when observed separately passes into inactivity. The 
same kind of notation is used in the description of the remaining transitions. 

Yet another typical unitary interaction in population models in the spontaneous birth of 
preys, E ^ A, but this is not considered here. All birth processes are due to those binary 
interactions that also occur in the system. At this point, it can be useful from a practical 
point of view to establish the probability v of considering rather a two-component transition 
than a single-component one. 

The first two-component interaction we will consider is the following: AB A EE. It 
may be specifically explained in terms of transactions among the qualified investors, but in 
a general sense, it conveys a form of agreement between two active investors in such a way 
that none of them convinces the other. This annihilation process is not usually considered in 
population models: it represents some sort of deadly defense in which predator and prey die 
after fighting. The usual result in predator-prey models after AB interactions is predation: 
AB A BB. In our case, this counts for the possibility that an active investor may perform 
a change in the evaluation of the market scenario (from bear to bull) due to predominance 
of S's. This may eventually lead to a market bubble. Once again it may become useful to 
consider that some fraction A of AB interactions conduces to annihilation, whereas 1 — A of 
them ends in predation. 

The third binary interaction is AE A/l, in which an agent that was not interested 
in operating in the market comes into activity in the A side. Under financial optics this 
imitative behaviour can lead to market panic and ultimately to a crash. Here lies our 
birth mechanism for preys that also incorporates in the model population pressure against 
unbounded prey growth. 

The complete state of the agent system at given instant t is fully settled by the number of 
investors belonging to species A and B, NA{t) and NB{t) respectively. Since these numbers 
will be stochastic magnitudes, we are interested in obtaining an expression for P(n, m, i), 
the probability of having n A's and m S's at time t: 

P{n, m, t) = Pr{NA{t) = n, Neit) = m}. 
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To this end we will consider the following five transition rates, the transition probabilities 
(per unit of time) between macroscopic states, based upon the above interactions: 



T{n — l,m — l\n,m 
T{n — 1, m\n, m 

T{n — 1, m + l\n,m 
T{n, m — l\n,m 
T[n + 1, m\n, m 
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n N — n — m 
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Note that there are also three "single-stepped" feasible transitions on each variable which 
are not listed here and, therefore, that are forbidden: T{n,m+l\n,m), T(n+ l,m + l\n,m) 
and T{n + l,m — l\n,m). 

Summing up, in essence we are assuming that: 

1. states A and B can spontaneously decay into inactivity, 

2. there is a basic non-trivial interaction (commonly the trade) that is not sensitive to 
the interchange of A and B roles, 

3. i?'s can convince A's only, and 

4. A's can convince E^s only, 

where the asymmetry in the two last items expresses the fact that bubbles and crashes 



in actual stock markets are different in shape 



2G 



21|. Since all agents are identical, the 



heterogeneity of our model relies on this asymmetry in the interactions. 



B. The master equation 

The Markov character of the model makes superfluous considering more sophisticated 
transition rates in the elaboration of the master equation (ME), the equation that defines 
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the time evolution of P{n,m,t): 
dP{n, m, t) 



dt 



= {aAA-lA){S^^ -l)[nP{n,m,t)] 
+ -fsiS^' -l)[mP{n,m,t)] 
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Here we have introduced the following increment /decrement operators 

Sf^f{n, m, t) = f{n ± 1, m, t), 



Sy^fin, m, t) = /(n, m ± 1, t), 



and five new parameters 

2t/c- (1 - u)p 

lA = , (1) 

IB = (2) 

o;aa = (3) 

Aa + (l-A)6 + c 
a^i? = 2z/ '-j;^ , (4) 

(1- A)&- Aa 

/i/AB = , (5j 

which encode all the relevant information of the model parameterisation. Let us stress that 
A and u were defined in order to clarify how the update mechanism can be approximately 
implemented, see Fig. [T], but they do not introduce further degrees of freedom to the prob- 
lem since they would disappear after a constant redefinition. This is the case if one uses 



Gillespie's exact algorithm [22] in the simulation of the system, as we have done. 

The relevance of the new parameters becomes noticeable soon afterwards. Suffice it to 
say for the moment that we will proceed as they were independent of the size of the system 
in what follows, because we will consider the expansion of the ME in terms of powers of A^. 
To this end let us define RA,B{t), 

RA,B{t) = lim EINaAWN. 

N^oo 



Otherwise 





+ St 






u~U{0,l) 



Pr{ AE ^ AA\ = cSt <— ^<choose 




Vr{A^E\=pSt 



Pi-{B-^£'} = 


qSt 




Pr{AB^£E} 


= aSt 



Figure 1: (Color online) Flux diagram for a discrete-time updating procedure of the state of the 
agents. 

and introduce two new stochastic processes, X{t) and Y{t), in such a way that 

NA{t) = N RA{t) + VNX{t), 
Nnit) = N Rsit) + ViVr(t), 

hold. X{t) and Y{t) are thus responsible for the fluctuations of A^a(^) and Nsit) around 
their mean values. It is expected that the strength of those fluctuations will diminish as the 
system reaches the thermodynamic limit, that is, when ^ 1. Note that this approach 
implies that, for any two given values of the species population, n and m, we will have that 

n = N RAit) + ViVx, 
m = N Reit) + ViVy, 

where x and y — as well as Ra{^) and Rsit) — are real magnitudes in spite that n and m 
were integers. In such a situation increment /decrement operators become partial differential 



operators 



23|, 



±1 



Finally note that P{n,m,t) must be replaced by Il{x,y,t), 

n(x, y, t)dxdy = Fi{x < X{t) ^ x + dx,y < Y{t) ^y + dy}, 
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P(n, m, t) = In r JL^\ JL^\ t ) dxdy, 



through 

what affects the time derivative term in the ME in the following way: 



dP 

dF 



1 ^^^/^TT^ 1 ^^S/^TT 



Vat dt " ViV dt " iV 

C. First-order stationary solutions 

The first-order approximation of the ME collects terms of order N~^^'^, ignores those of 
0{N~^), and leads to a set of coupled Volterra equations for RAit) and Rsit), 

dRA 



dt 
dRs 



[lA - ^aaRa - aABRs] Ra, (6) 
[/3abRa-1b]Rb. (7) 



dt 

Let us analyse the factors appearing in these equations. 'Ja as defined in Eq. ([T]) represents 
a trade-off between a positive term that comes from the imitation influence and a negative 
term that measures the death rate of preys. If positive, it would correspond to a effective 
birth rate of preys in Eq. ([6]). Recall however that in this system preys suffer of population 
pressure instigated by the imitation interaction that constrains the preys' growth, see the 
definition of oaa in ®- The term with the oab factor counts for the reduction in the prey 
number due to all binary operations, not only predation, Eq. (jl]). The f3AB term appearing 
in Eq. ([7]) is a consequence of the imbalance between predation and annihilation alternatives, 
as it can be observed in ([5]), whereas •jb measures exclusively the death rate of predators, 
expression ([2]). Summing up, there are two parameters, 'Ja and Pab, with no definite sign, 
whereas 7^, oaa and oab are positive constants ab initio. 
Eqs. (E]) and ([7]) present three stationary solutions for which 

dRA dRe 



dt dt 



0. 



The first solution is the trivial one, Ra = Rb = 0- It represents the death of the market 
due to a complete lack of activity. This is a feasible scenario that threatens any real market. 
For instance, investors may loose interest in any given commodity that becomes useless or 
exhausted. The stability analysis of this fixed point determines that it will be a saddle point 
if jA > 0, otherwise it would turn stable. The analysis of the second stationary solution, 
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Ra = 7a/«aa = M/N < 1 — note that 7^ < aAA by construction, c/. expressions ([T]) 
and dS]) — , and i?^ = 0, leads to the constraint 

if one wants to avoid conferring stabihty to this fixed point as welL ^ In conclusion, all the 
parameters defined in dl])-® must be positive-definite. 

We must point out that the presence of those unstable equilibrium solutions is not a flaw 
but a merit of the model, as is the fact that the remaining stationary solution 

D DO — 1b 

J^A — J^A 



Rb — R°B = 



Pab 

IaPab - iB^AA 



o^abPab 

is always present, accessible and corresponds to a stable fixed point. 

Regarding the occurrence of the fixed point, it is evident that R\ > 0, and Eq. (IHl) leads 
to R°j^ < M/N < 1. The same equation determines that R% > 0. Also R°^ < 1, as it can be 
proven as follows: 

DO 1 ("AB - 1a)Pab + IBOIAA ^ 

^B = ^ n < J-) 

OiABPAB 

because trivially aAB > 1a, cf. Eqs. ([T]) and @. We can also show that i?^ + R°^ < 1, 
DO , DO 1 ("AS - 1a)Pab + (aAA - oiab)ib 

^A + ^B - -L n 

aABP. 



< 1 



AB 

iB {aAB - lA)aAA + (qaa - (^ab)1a 

lA OiabPab 



^ ^ _ lB{aAA-lA) ^ ^ 

IaPab 

because aAA > 1a as we have just pointed out above. 

The analysis of the stability of this fixed point leads to the conclusion that the point is 
stable, and that the transient term will exhibit oscillations when uoq G M^, 



c^o = \l uabI^abRaRI - \ {c^aaRa?, (9) 



^ Note how this situation may represent the worst market crash imaginable, where every active investor 
is in the bear side, what would conclude in a liquidity crisis like the present one. In this sense one may 
consider that M is related to the finiteness of the total amount of shares. 



11 



"^^<2, 1 + ^-2. 



which is true whenever 

When the system shows transient oscillations, there is a single characteristic time scale for 
the decay rate, 

and for t ^ tq the system would reach the stable solution. This assertion is no longer true 
when a second decay rate appears, since 

To"' = Tq-^ - aABPABR^RB < ro^- (11) 

In such a situation the steady state is reached when t^^ -C Tq'^ — Tq"^. Therefore we may 
define t^, t^^ = Tq^ — 3?[Tq"^], and the steady state is always achieved for t ^ tg- 

After the transient regime, and whenever is finite, we will await that the time evolution 
of prey and predator densities, NA{t)/N and NB{t)/N, makes them attain their limit values 
and i?^, and exhibit some fluctuating activity afterwards. Since the characteristic size 
of the fluctuations is of order N~^^'^, a naive analysis could lead to the conclusion that if we 
have, let us say, 1000 interacting agents the error in neglecting the remaining terms in the 
ME should be around 3.2%. In Fig. [2] we can find the outcome of a realisation of the model 
with = 1000 — the complete set of parameter specifications is listed below in Section [XTTl 
The example shows that in such a system fluctuations may be larger than expected, and 



further corrections to the first-order equations must be taken into account 10, l24|. 

A final word on the A^ dependency of the above expressions before exploring the forth- 
coming terms in the ME. The order-by-order analysis under progress relies on the fact that 
the parameters defined in Eqs. ([I])-® are independent of A^. Note however that the expres- 
sions for R°j^, R°Q and M/N are insensible to this need. It only affects those constants where 
time is involved, like tq or c^o- 
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Figure 2: (Color online) Time evolution of preys and predators densities (solid red line) for an 
exact realisation of our interacting agent model with N = 1000. The dashed black line depicts the 
first-order approach to the problem, whereas in green is shown the stationary solution. We can see 
how fluctuations in both populations are larger than expected. 

D. Beyond the first-order equations 



When one gathers the terms of order ^ in the ME expansion, a Fokker-Planck equation 
for Il{x,y,t) emerges: 

dtU = [-'yA + '^aAARA + aABRB]dx{xU) 

+ aABRAdxiyTl) - pABRsdyixU) + [jB - PabRa] dyiyU) 
+ -Y [-^^ + "^^(2 -Ra- 2Rb) + aABRs] d^.U 



R 



+ ^ [7b + ("AS - aAA)RA\ d''yyU - PABRARBd'iyU. 

If we concentrate our analysis of the previous equation for times large enough to let RA{t) 
and Rsit) reach their steady state values, R°^ and R^, the expression simplifies considerably: 

dtU = fx^AixU) + fJ^xydxiy'n) - l^yxdyixU) 
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positive-definite quantities. ^ Therefore we have a linear multivariate Fokker-Planck equation 
for the joint probability density of X{t) and Y{t), whose solution can be systematically 
obtained after some (or maybe plenty of) algebra [23]. An alternative approach is based in 
the following set of coupled (Ito) stochastic differential equations: 

dX = -fixxXdt - fXxyYdt + axdWi, (12) 
dY = fXyxXdt - paydWi + Oy^l- p^dW2, (13) 

where W\ and W2 are two independent Wiener processes. 



E. The magnifying effect 

In order to explore the reason for the abnormal magnitude of fluctuations we should com- 
pare Xii) and Y{f) with R\ and R°^ respectively. A quick analysis reveals that mean values 
are not useful in this venture because limf_,oo IE[X(t)] = limt_^oo E[F(t)] = — remember 
that (|T2l) and (|T3l) are valid for t ^ to- We concentrate in variances and co- variances instead. 
In Appendix |A]we can find how the stationary values of E[X2(t)], E[y2(t)], and E[X(t)F(t)] 
follow: 

a.(o) = iimE[x^(t)] ^y^'^'^^^y^ 



'^Pxxpyx 

Pyx'^x ^ (pxx ~^ pxyPyx) Cj, ~ '^PPxxPyx'^x'^y 



CyyiO) = hm E[Y\t)] 

t— >oo ^PxxPxyPyx 

2 

Cxy{Q) = lim E[X(t)r(t)] - 



yx 



Note in particular that j3AB + chab — o^aa = 4i^(l — X)b/N > 0. 
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and fluctuations will be typically large if these quantities are much more bigger than (RX)^, 
{R°bY s-^d ^°A^°By respectively. If we deflne the magnifying factors fl^x, ^yy and fl^y as the 
corresponding quotient of these magnitudes, e.g. 

E[X{t)Y{t)] _ a,(0) 



n 



xy 



lim 

t— >oo 



Do Do 



fluctuations will be relevant when Q ~ A^, because then one can overcome the N'^l"^ dumping 
factor of the second-order corrections. The analysis of the possible values that can take 

is difficult because the inner relationships that ^^x^ fJ'xy, fJ'yx, o'x, Uy and p present. In fact 
the difficulty is inherited from 7^, 7^, c^aa, chab and Pab, which are neither bounded nor 
independent. Then, it is useful to introduce the following (final) re-par ameterisation: 

1 2 



OlAA = 

OtAB = 

Pab = 

lA = 

IB = 



1 2 



1 



1 + 

X 

l-r)2 



X 



2 

To 



V To 

where the four new variables ^, V and ^ are in the (0, 1) range, and can be arbitrarily 
settled. With the proposed parameterisation all the constraints that affect the old param- 
eters (the pure algebraic ones, as well as those coming from stability considerations) are 
identically satisfied, ^ and tq carries the characteristic time scale of the interactions at the 
microscopic level. The magnifying factors in the new parameters read 

1-x 1 



;i -rye)- 



X^ 



2 ^ XV 



yy 



1 



X, 



+ 



X 

+ e 



1 - r^e) - 1 
X 



1+e 



1—7] 



X)rfe 



2 e (1 

+ e 



1 



xy 



2 e (I-X)^e' 

and the stationary first-order solutions are i?^ = x and = (1 — x)v^- 



The new parameterisation proves incidentally that p < 1/2. 
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The first point to be noted is that no tq appears in any of the last expressions. So, the 
magnification effect does not depend on the characteristic time scale of the correlations. 
The second aspect of importance is that, for fixed values of ^ aiid rj, the magnifying 
factors become unboundedly large as ^ —^0, and this parameter does not contribute to the 
value of R°j^ and R°q. Then, magnification can be achieved for any (regular) value of the 
species stationary densities. Another favourable scenario is i?^ and R"^ — * 0: check for 
instance how for x — > 0, Q^x ^ oo. This implies that the phenomenon is relevant in sparse 
systems as well, in spite of the fact that in such cases may eventually be very large. Note 
finally that magnification is not connected with the presence of oscillations of any peculiar 
frequency. On the one hand, the condition that determines that Tq~^ replaces ujq is 



4(1 -x)(l - >;)£' 

and as we have shown above ^ — > always leads to magnification. This is reasonable since 
for a fixed tq, Tq increases the microscopic correlation range — see Appendix |Al On the 
other hand, for a fixed value of R°^, R°q and ^, ujq embraces the whole positive real axis as r] 
varies. Therefore, in principle, one can consider models with either a large value for cuq, and 



reproduce the typical bid- ask bounce in a liquid market [25|], or a smaller one, and capture 
some seasonal character present in the market evolution, like in electricity markets (26 1. 
Because, as it is shown in detail in Appendix El and it can be observed in Fig. [21 the 
oscillatory behaviour is also present in the second-order terms. 

Let us see magnification in a practical example. For clarity reasons we will condense the 



three magnifying factors defined above in a single plot. To this end we define VL^ 

^zz — ^xx ~l~ ^yy "^^xy lim E 



t—*oo 



m_xity\- 

R% R\ ) 



a relevant quantity, as we will see below. Further, we will fix R\ = x, and change e in a 
way that keeps R°^ constant, what leaves t] and ^ as the only free parameters. In particular, 
we have set R°^ = R°q = 0.2, since we are interested in models in which the two sides are 
equally (no side is prioritised) and macroscopically populated. In Fig. [3] we observe some 
contour lines that represent configurations with the same amplification level, and how these 
lines cross the borderline that delimits those configurations with and without oscillating 
properties. Thus, for instance, we have marked with a small circle the location of the 
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following parameter set: x — 0-2, e = 0.625, 77 = 0.4, ^ = 0.2. With this configuration 
M — 0.7 N, and the amphfication factor is about one hundred. 




Figure 3: Contour plot of the magnifying factor Qzz- The magnifying factor values, to be compared 
with N, are 100 (dotted line), 250 (dashed hne) and 1000 (dot-dashed line). The solid line is the 
borderline between the zones with oscillating behaviour (wq) and without it (To). 



III. PRICE DYNAMICS IN A EXCESS DEMAND MODEL 

As we have stated above, the formula that will determine the price dynamics must depend 
on the nature of the species. Then it is time to identify A and B states and define how 
the evolution of the population of agents in each category translates into prices changes. 
Within this first model we will generically consider that A^a(^) represents the offer and 
Nsit) represents the demand in a certain financial market, and assume that in such a 
market excess return reacts linearly to excess demand: a classical and ubiquitous point of 
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view in the economic literature, e.g. {4, 11, 14, Q, I^]. Excess return measures the 
logarithmic earnings of the stock beyond the risk-free interest rate r, R{t) = In [S'(t)e~''*], 
and excess demand is the difference between Nsit) and NA{t). Therefore we will have 

dR{t) = ^{Nb - NA)dt 



where S measures the sensitivity of prices to excess demand. This first pricing model is 
a good testing ground since the agent model will be responsible of any observed market 
property: we are simply integrating the differences in population. 

Let us consider the following paradigmatic example with = 1000 investors. We have 
set To = 10 minutes, so it is of the same order of magnitude as a typical correlation length 



found in actual financial data [28j. Beyond this, the rest of values were not based on actual 
market observations. In fact we have set x = 0-2, V = 0.4 and ^ = 0.2, like in the example 
we emphasised in the previous Section, but slightly increased the value of e, e = 0.643. This 
was intended to get R% > R% R% ^ 0.206, whereas R\ = 0.2. Note that R% - R'X > 
characterises a growing economy in which wealth is injected into the market. This term is 
also responsible for a long-run exponential growth. 

A possible realisation of the dynamics of the species population was previously introduced 
in Fig. [21 and in Fig. H] we find the corresponding evolution of the stock price when S = 10~^ 
min~^. Here we sampled the complete data series to consider closing prices only, a usual 
practice in technical analysis. Moreover, in the confection of Fig. H] and hereafter we are 
assuming that a trading day lasts 480 minutes, and that there are 250 trading days in 
a year. We observe in Fig. H] the appearing of typical market charts: upward trends — an 
increasing succession of minima — , downward trends — a decreasing succession of maxima — , 
and sideways trends — a bouncing movement between two price levels. 

In Fig. E] we present the outcome of a statistical analysis performed with the stationary 
data set of fixed-time returns R{t] t) = R{t + r) — R{t), t > 100 minutes. We check that for 
r ~ To correlations are important, Gaussian limit is not attained and skewness is observed. 



like in actual markets 



28 



29 



30 



3l| . This phenomenon is even more noticeable when 
the standard deviation of fixed-time returns, a measure of the volatility of the market, is 
analysed. Fig. [HI Since X{t) and Y(t) are anti-correlated, and the return change is sensible 
to the difference of those magnitudes, we expect that volatility grows faster for small time 
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Figure 4: (Color online) Time evolution of the daily closing value of (discounted) stock prices. We 
can see how the market suffers a market bubble (an upward trend followed by a downward trend) 
lasting 5 years and ending at the beginning of year 25. After that we find what is called a sideways 
trend, i.e. no trend at all, from the mid of year 25 to the mid of year 27, followed by a new upward 
trend with corrective movements in the middle. The inset shows the exponential growth in the 
long run. 



scales, and reaches the diffusive regime for t > tq. Traits of abnormal 



both sub- and super-) 
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33|. 



diffusion have been reported to be present in real markets as well 

In the confection of the previous plot we have used the complete set of returns available 
for each time scale r, by assuming the statistical equivalence of every sample R{t; t) as a 
function of t. Moreover, the results above seem to indicate that, for r ^ tq, the samples 
R{t] t) and R{t] t + t) ought to be also (almost) independent one from the other. So, if we 
compute the realised n-r volatility, Ki(r; k): 



Vn{r\ k) 



1 r 1 

i?(r;(A;-m)r) - -i?(^nr;(A;-n)r 

\ m=l L 



we should expect the outcome to be uniform in fc, as well as an absence of correlation between 
Vn{T\ k) and V^(r; A; + n). In order to check whether this assumption is true we have chosen 
r = 1 day, and = 20 trading sessions, as a proxy for the one-month realised volatility, a 
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Rescaled returns 

Figure 5: (Color online) Fixed-horizon return behaviour. We can see how probability density 
functions at small time scales slightly but noticeably depart from Gaussian behaviour by exhibiting 
negative skew: the negative tail is fatter than the positive tail. Returns were divided by their 
sampling standard deviations to make them commensurable. 



typical choice among practitioners. ^ The results were also annualised, which means here 



that they were increased by a factor a/12.5, and only k ^ 21 are considered — we ignore the 
whole first day of simulation. The outcome, as it can be observed in the inset of Fig. [TJ is that 



the market alternates long periods where the vo 



a phenomenon known as volatility clustering 3l|| • The presence clustering in the volatility is 



atility is large, with periods of relative calm, 



a well-documented feature of real markets that is usually explained in terms of the existence 
of volatility self-correlation. This correlation, as opposed to the return-to-return correlation. 



is long-ranged 



34( 1 — confront time scales in Figs. [7] and [12] below. 



In order to offer a plausible origin for this larger time scale, we have composed Fig. [HI 
There we present, in a phase diagram, the mean recurrence time: For each possible state 
of the system we have recorded all the visiting times, and performed a sample mean with 
the inter-event times. Therefore, at least two visits to a given state are needed in order 



^ For instance, one-month volatilities are on the basis of VIX, the Chicago Board Options Exchange pio- 
neering volatility index for the Standard & Poor's 500 index. 
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Time (minutes) 

Figure 6: (Color online) Volatility growth. In this figure we can see how the volatility growth 
presents two well different regimes. For r < tq the standard deviation of fixed-time returns shows 
(super-diffusive) linear growth whereas for r > tq it scales as -y/^) like in a diffusive process. 

to attach a non-zero value to that point. Once again we have disregarded the data within 
the first day. As it can be observed in this figure, the mean time grows in an exponential 
fashion when we depart from the stable fixed-point values, Na = 200 and Nb = 206. Since 
the scale is logarithmic some lower bound must be chosen, and we have decided to remove 
those data points with a mean recurrence time smaller than 1 minute. This retains in the 
plot almost all non-zero values: ^ the lowest recurrence time near the core is attained at 
Na = 207 and Nb = 203, yielding a value of 20.17 minutes. As we can see in Fig. [T^ again, 
this magnitude coincides with the time scale for which one-minute returns exhibit stronger 
anti-persistence. The slow decay in the volatility self-correlation may have thus its origins 
in those long periods needed by the system to return to the most outer zone, from where 
the largest absolute returns come: the light (green) ring marks a recurrence time of about 
60 days, the time scale for which the volatility self- correlation is more intense — see Fig. [71 
Therefore, the agent model is capable of conciliating short-range correlations in the mi- 

^ Indeed, with this practice only few peripheral points were ignored: states that were visited several times 
in a rapid succession before the system leaved that zone and never returned to it. 
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Figure 7: (Color online) Volatility clustering. In this figure we can see how the one-month annu- 
alised volatility presents clustering (inset) and long memory. The increment of the correlation for 
times below 20 days is due to data overlapping. 

croscopic level with long-range correlations in the macroscopic level, what may be linked 
with the presence of large business cycles in the financial data 351]. 

Finally, another stylised fact that is commonly associated with clusterin^and l ong -range 



memory in the volatility is the so-called leverage or Fisher-Black effect |31l . l36l . l37l |. This 
phenomenon is generically characterised by a negative relationship between returns and 
volatilities. In our case, this effect can be barely observed when the cross-correlation between 
20-day returns and volatilities are depicted — see Fig. [9l We must point out, moreover, that 
there are features shown in empirical studies related to this effect that are not detected in our 
example. For instance, from Fig. one cannot sustain the presence of a noticeable temporal 
asymmetry in the correlation, as expected [37]. However, we can explain this departure from 
what is observed in actual markets on the basis of the usual interpretation of the leverage 
effect: the market digests with nervous the loses, and with confidence the rises. And we 
must remember at this point that the price information is not fed back into the species, so 
it is not possible such a reaction here. Therefore, the slight anti-correlation present in Fig. [9] 
could be a side effect of the volatility self-correlation, but just some spurious result as well. 
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Figure 8: (Color online) Phase diagram of the state of the system. We show the mean recurrence 
time, the mean time passed between two consecutive visits to the same state. 

IV. PRICE DYNAMICS IN A LIQUIDITY MODEL 

Another possible identification of the two species comes from considering that an agent 
in the A state is a liquidity provider, whereas an agent in the B state represents a liquidity 



taker 
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381]. Liquidity providers introduce hmit orders into the market whereas hquid- 



ity takers operate through market orders. Limit orders are orders with a hmit price that 
represents the minimum (respectively maximum) price the investor is accepting for selling 
(respectively buying) a given number of shares, the volume of the order. The limit order is 
placed in the so-called limit order book, which is visible to rest of qualified investors and it 
remains there until one of the two following major events takes place: someone accepts the 
ask (respectively bid) price and the transaction is completed, or the investor removes the 
order from the book. Market orders are orders that automatically match the best opposite 
limit order out of the limit order book. This is just the way in which noise traders are as- 
sumed to participate in the market, but unlike them, liquidity takers survey the limit order 
book and look forward to a convenient trading opportunity, like predators do 39|. Note 
however that this results in transaction {AB EE), not in predation {AB BB). This 
and the rest of interactions determine how qualified agents migrate from one set to the other 
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Figure 9: (Color online) Leverage effect. In this figure we can see how 20-session returns and 
volatilities are negatively correlated, a phenomenon known as the leverage effect in the literature. 



or remain inactive 40l. l41|. 

This interpretation will be more relevant in liquid markets, where changes in offer or 
demand have a small impact on prices. In such a situation a relevant magnitude in the 
price formation mechanism is the spread, the difference between the lowest ask price and 
the higher bid price. Then, the pricing formula must connect the spread with the relative 
populations of liquidity providers and liquidity takers, but the issue is not so unambiguous 
as in the previous case. We have decided to use a pricing expression inspired by [42, l43 |. 
Consider these general guidelines: 

1. The bigger the number of limit orders are in an order book, the lower the spread will 
be, and therefore, the lower the price change will be. 

2. If the number of liquidity takers is small with respect to the number of liquidity 
providers, the price should tend to exhibit the typical bid-ask bounce pattern. 

3. If the number of liquidity takers is large with respect to the number of liquidity 
providers, the most likely is that price shows a trend. 

A feasible candidate that incorporates the above properties is the following discrete-time 
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updating formula: 



R{t + At) = R{t) + EG {R{t) - R{t - At)) 



Nsit) 
NAit) 



C 



At, 



where At is the time between two consecutive changes in the agents' state, 0(-) is the 
Heaviside step function and, for a matter of model simplicity, we will assume that ( = 
R°q/ R\. This liquidity model will share some traits with the previous excess demand model, 
since for large values of N and t we have 

-Y{t) X{t)- 



R{t + At) ~ R{t) + {R{t) - R{t - At)) 



N 



R1 



R\ 



At, 



and the two formulas become very similar when R°^ ^ i?^. The main distinguishing trait is 
the presence of the factor with the Heaviside step function, that may distort the evolution 
that the agent model dictates and introduce new properties. 

In Fig. [TUlwe see how a sample time series mimics again a typical stock market evolution. 
In the confection of this plot we have kept the same parameters as in the previous market 
model with just one exception, the sensibility was set to S = 0.05 min~^, with the aim of 
recovering a similar growth in the long run. However, we see how the market becomes much 
more volatile than it was in the previous case. 

The impact of the new pricing dynamics in the probability density function of returns is 
also noticeable. In Fig. [TT] one observes how one may find in practice changes amounting 



tens of standard deviations, like in actual markets |28l. I29l| . Finally, the correlation length is 
affected but time scale tq is preserved: in Fig. [12] we show a comparison between the one- 
minute-return correlation curves for both cases. On the contrary, the pricing mechanism in 
the liquidity model wipes completely the negative correlation that one can relate with ujq: 
the characteristic time scale of the oscillations is about 44.4 minutes with our parameter 
selection. 



V. CONCLUSIONS 

Along this article we have introduced a dynamical model that describes ultimately the 
evolution of financial prices. The main ingredient of the model is a finite set of identical 
interacting agents that at every moment can be accommodated into one of three excluding 
categories. The agents represent those traders whose operations may noticeably impact on 
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Figure 10: (Color online) Time evolution of the daily closing value of (discounted) stock prices for 
the alternative dynamics based in the liquidity model. Among the several trends that are present, 
it is remarkable the sudden large drops which resemble market crashes. 

the market, and the three available states characterise in a broad sense the investor possible 
attitude. 

Active agents can spontaneously adopt a neutral attitude, but any other change is the 
outcome of an agent-to-agent interaction: agents may agree, or one agent can convince the 
other following hierarchical relationships. These simple rules encode a system in which the 
number of active agents may strongly fluctuate, thus overcoming the second-order nature of 
the effect. We have looked for the conditions that promote such amplification, and concluded 
that it does not depend on the time scale of the interactions, can be obtained for any choice 
of the first-order stationary densities — even though it is more relevant in sparse systems — , 
and it is not the result of a resonance. 

Once we have analysed the dynamics of the agent instantaneous properties, we have 
moved into the pricing problem. We have considered two different identification of the 
categories and, in each case, a suitable pricing expression is set. We have simulated the time 
evolution of the asset price for representative values of the involved parameters. We have 
shown how sample realisations reproduce several stylised facts reported in actual financial 
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Figure 11: (Color online) Fixed-horizon return behaviour. We can see how probability density 
functions present fat tails at intra-day time scales, and how the Gaussian behaviour is not fully 
recovered even in the case of daily returns. 



data sets: the price evolution displays upward, downward and sideways trends; probability 
density functions of small time-scale returns present fat tails and skewness; volatility behaves 
accordingly in a non-diffusive way within the same time horizon, and presents clustering in 
a larger time scale; and traces of some leverage effect can be found. 

In a forthcoming work we are planning to explore how the properties shown by the agent 
model depend on the assumptions made, to refine its connections with actual financial 
systems, and to consider further alternative interpretations that may be relevant in market 
dynamics not considered here. 
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Figure 12: (Color online) Linear correlation of one-minute returns for the two pricing models. The 
basic time scale tq = 10 minutes can be observed in both cases. 

Appendix A: CORRELATION FUNCTIONS 

The cross-correlation theorem states that, for any couple of two random variables X(t) 
and Y{t) we can compute its stationary auto-correlation function though: 



X{uj) stands for the Fourier transform of X{t), and so forth. When Y{t) coincides with 
X{t), Pxx{uj) is termed the power spectral density function of X{t), and cross-correlation 
theorem is known as the Wiener- Khinchin theorem. 
In our case equations (|T2l) and (|T3l) lead to 




where 
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f^xxf^xy^y ~l~ f^xyf^yxP^x^y 
^^/^L + (^^ - PxyPyxf 
i [Pyxf^l + Pxy<yl - PxxP<yx<yy) ^ " PCr^O"? 



^'^P'lx + (^^ - PxyPyxf 

Note that every function has the same basic structure, namely 



Piuj) 



^^/^L + i^^- PxyPyxf 

therefore we can compute the auto- and cross-correlation functions at once, 



C(r) 



^1 + n^PxyPyx ^ ^1 /^3/^x,/i,. sin(u;ok|) 



PxxPxyPyx 



'^PxyPyx^O 



+ 



K2 



Pxx^O 



■ sm[uJoTj 



■\r\/To 



where uq G M+ coincides with that defined in 1^, and tq was introduced in (ITUl) . If Tq G 
see Eq. ffTTj) . we have instead 



C(r) 



/^l + f^SpxyPyx 
PxxPxyPyx 



cosh (r /To) 



fiSpxyPyx 
'^PxyPyx 



Tosinh(|r|/To 



H Tosinh(r/To) 

Pxx 



■\t\/to 



and two correlation time scales appear. Finally note that the corresponding variances and 
co-variances are obtained after setting r = 0. 
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